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It has recently been established that hadronisation corrections to QCD jets vary as 1/R, 
at small R, for jets of radius R. Here we demonstrate, using jets in the k t algorithm, 
that the magnitude of these 1 / R corrections are unambiguously linked to the magnitude 
of 1/Q corrections to commonly studied event shapes in e + e~ annihilation. 

1 Introduction 

An understanding of QCD jets and their properties will be integral to the success of the LHC 
physics program. In particular one of the most important issues in current jet studies is the 
question of jet energy scale. The shift induced in jet energy by effects such as perturbative 
radiation and non-perturbative effects like hadronisation and the underlying event would 
contribute to a smearing of, for instance, mass peaks that may be a signal for new physics. 
Thus in order to choose optimal jet definitions that minimise such smearing one would need 
to know the dependence of the above effects on the experimental parameters like jet radius. 
Moreover, even in pure QCD studies such as the extraction of parton distribution functions 
(pdfs) and the strong coupling a s from jet observables like the inclusive jet cross-sections, 
a knowledge of the non-perturbative contribution is important to supplement perturbative 
calculations. A relatively small shift in the transverse momentum p t of a jet, induced by 
hadronisation, can result in a significant change in the inclusive jet spectrum since we are 
dealing with a quantity that has a steeply falling p t distribution. 

While it is traditional to study the hadronisation contribution via Monte Carlo models 
such as those in HERWIG and PYTHIA, it turns out that in cases like the jet energy there 
is additionally valuable analytical insight available [2J. Analytical models based on renor- 
malons [3] have in the past been met with great success in the description of LEP and HERA 
event-shape variables [4], but have not really been utilised outside that context. In Ref. [2] 
one such model (due to Dokshitzer and Webber [5]) was used to estimate hadronisation cor- 
rections to jet transverse momentum pt- The result found there was striking: hadronisation 
effects have a singular 1/R dependence on the jet radius R, at small R. This is in complete 
contrast to the contribution from the underlying event which varies as R 2 . The knowledge 
of the R dependence of non-pcrturbativc effects in conjunction with the In R behaviour in- 
volved in perturbative estimates can then be used to arrive at conclusions about the optimal 
values of R to be used in diverse studies involving jets, as exemplified in Ref. [2J. 

While the computations of Ref. [2J indicate the dependence of hadronisation corrections 
and the underlying event on R, there remains the question of the overall magnitude of these 
effects. While for the underlying event one is reliant solely on Monte Carlo event generators 
to obtain the overall magnitude, for the hadronisation correction a tentative link was made 
in Ref. [2] between the magnitude of the 1 /R correction and that of 1/Q corrections to LEP 
and HERA event shapes such as the thrust distribution (see [4] for a review). In order to 
definitely link the magnitude of jet hadronisation to that of event-shape power corrections 
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one needs to carry out a calculation at the two-loop level rather than the simple one-loop 
estimate reported in [2]. The calculation for jets defined in the k t algorithm [SI [7] is reported 
here while the corresponding result for jets in the anti-A:* algorithm is already known [8]. 
Work on the other jet algorithms is currently in progress. 



2 The single-gluon result 

In the Dokshitzer- Webber model non-perturbative hadronisation corrections are associated 
to the emission of a soft gluon with transverse momentum k t ~ Aqcd- For the jet pt case 
we work out the change in transverse momentum Spt induced by the emission of such a 
gluon and combined with the gluon emission probability (as given by perturbative QCD) 
this yields the average shift in pt induced by hadronisation: 



where kt, T] and <fi respectively denote the transverse momentum, rapidity and azimuth of 
the emitted gluon with respect to the emitting hard parton (jet) and Cj is the colour factor 
associated to emission from the given hard parton. The reason we are able to single out a 
given hard parton initiating a high-p t jet in any hard process is essentially since the leading 
1/R result stems from emission collinear to the triggered jet. It is thus possible to talk in 
terms of the single-jet limit ignoring the rest of the details of the hard process. 

The only non-perturbative ingredient that is involved above is the value of a s (kt) at 
scales around or below Aqcd- If one makes the assumption of a universal infrared- finite 
coupling, which replaces the perturbative coupling that has an unphysical divergence at 
Aqcd, then one arrives at a prediction for the hadronisation correction (Sp t ) h . Using the 
fact [2J that the 5pt(k) is essentially the energy of the gluon emitted outside the jet we 
perform the integral over rapidity in Eq. |T]) to obtain the leading 1/R behaviour. The 
result is of the form cA/R, where c is a number obtained from the rapidity integral and 
A is the moment of the coupling a s {k t ) over the infrared region (we refer the reader to 
Refs. [SlUn] for the precise details). Since the same coupling moment enters the predictions 
for event-shape variables we can take its value from data on event shapes and hence obtain 
a numerical prediction for the leading 1/R hadronisation correction to jet pt- This was the 
method adopted in Ref. [2j. 

Here we point out a limitation of the above approach which is that while we have 
written down and used a running coupling a s (k t ), this quantity only emerges when one 
considers not just the emission of a single gluon but in fact gluon decay as well. To be 
precise an inclusive integration over gluon decay products is responsible for building up the 
quantity a s (k t ). Unfortunately, as is known for event-shape variables, our observable is 
sensitive to the precise details of gluon branching and hence one is not free to carry out such 
an inclusive integration. One must therefore return to the details of the gluon branching 
and identify the correction to the above inclusive approximation. The analysis at this level 
has already been carried out for event-shape variables O [101 H2J H3] and below we report 
on it for the jet pt case. 
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3 Non-perturbative effects and gluon decay 



Now we consider the situation where the emitted gluon with k t ~ Aqcd is allowed to decay 
and at the same accuracy account for virtual corrections to single gluon emission, as depicted 
in Fig. [TJ 

At this two-loop level the change in p t (for a quark jet) can be expressed as: 



(Spt 



7T 



d 2 k t da 
wk 2 a 



{a s (0) + 47r X (fc t 2 )} 6pt(k) + 4C F 



2 m 2 
dT 2 -^-6p t (ki,k2), 




where a is a Sudakov variable, a s (0) is an 
ill-defined quantity which will cancel away 
subsequently, x represents the virtual cor- 
rection to gluon emission, dT2 is the gluon 
decay phase-space and M 2 is the decay ma- 
trix element [U QU]- We also denote by 
Spt{ki, k 2 ) the change inp t due to correlated 
two-parton emission while 5p t (k) is the cor- 
responding singlc-gluon quantity. To cor- 
rectly account for gluon branching one thus 
has to perform the above calculation, the 
details of which are reported in Ref. Q3] • 

The analogous two-loop analysis for 
event-shape variables [9l QUI 021 H3] revealed an initially surprising result - the two-loop 
correction simply provided a universal rescaling factor to the one-gluon result, which be- 
came known as the Milan factor. Its value for nj = 3 (which is the number of flavours 
excited in the relevant soft region) was found to be M. = 1.49. Thus the ratio of corrections 
to two event shapes vi and v 2 was merely the ratio of the one-loop coefficients computed 
previously [5]: 



Figure 1: Gluon decay and one- loop correc- 
tions to singlc-gluon emission. 
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where Ji;^' 1 denotes the non-perturbative single-gluon correction for v\ computed as dis- 
cussed in the preceding section and likewise for v 2 . This remarkable result was understood 
to arise as a consequence of the fact that all the variables considered could be expressed as 
linear sums over the transverse momenta ku of emissions, v = ^2 i k t iCi, where the Cj are 
rapidity-dependent coefficients [TU] . 

In the case of jet pt this linear dependence is ruined by the non-trivial action of the 
jet algorithm in all cases except the case of jets defined in the anti-fct algorithm [8]. The 
contribution to the jet Spt of a given emission is found to be of the form /^e^E^ut, where 
•q denotes the rapidity with respect to the emitting hard jet and 2 out denotes the condition 
that the emission ends up outside the jet after the application of the jet algorithm. It 
should be immediately clear from this that in most current sensible jet algorithms (both 
of sequential recombination and cone type) the condition S out is non-trivial and introduces 
non-linearity in k t . For instance in the k t algorithm we can consider the situation in Fig. [2J 
where although one may have a soft parton separated by more than a certain distance R in 
rapidity and azimuth (denoted by the red gluon line) from a given hard parton, it may be 
clustered to another soft parton (denoted by the black gluon line) and hence swept into the 
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final jet. This clustering depends on the k t of a soft parton relative to the other partons and 
hence the condition S out derived in |14j contains dependence on the k t of the soft partons, 
spoiling the simple linear dependence needed for universality. 

An exception to the above situation is to 
be found in the anti-^t algorithm for which 
the condition S out = 0(?y 2 +</> 2 — R 2 ) ensures 
that a given parton is outside the jet if its 
angular (77, 0) separation from the hard jet 
is more than R. The linear dependence on 
kt is maintained and the Milan factor M. = 
1.49 is computed as for event shapes. 

For the k t algorithm we have carried 
out an equivalent calculation for the lead- 
ing l/R hadronisation correction (at small 
R) with the more complicated S out func- 
tion involved there and we found the result 

Mk t = 1-01 (rif = 3). Thus while at the level of the one-gluon studies of Ref. [2] the kt 
and anti-fct algorithms received identical hadronisation corrections, a detailed analysis at 
the two-loop level breaks this equality. One finds that the ratio of hadronisation corrections 
is then: 




Figure 2: The role of clustering in determining 
whether a soft gluon ends up within or outside 
a hard jet. 
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Thus one expects somewhat smaller hadronisation corrections for the k t algorithm as 
compared to those for the anti-fc t algorithm which is also borne out by the Monte Carlo 
studies with HERWIG and PYTHIA reported in Ref. [2J. We remind the reader that these 
conclusions apply only to the l/R hadronisation corrections that would be dominant at 
small R and we neglect finite R corrections which need to be considered alongside the 
underlying event contribution which also has a regular R dependence ~ R 2 . 



4 Conclusions 

We have reported on a study of hadronisation corrections to jet pt or energy scale based on 
two-loop extensions of the one-gluon estimates reported in Ref. [2], with jets defined in the 
kt algorithm. Studies for other jet algorithms (SISCone [15] and Cambridge/ Aachen [16]) 
are in progress. 
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